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SUMMARY 


Theoretical analyses are made of the stress distributions in a 
sheet of parallel filaments which carry normal loads and are imbedded 
in a matrix which carries only shear. In all cases, uniform loading at 
infinity is assumed and small- deflection elasticity theory is used. 
Static and dynamic stress- concentration factors due to one or more fila- 
ments being broken are determined. Particular attention is paid the 
dynamic overshoot resulting when the filaments are suddenly broken. The 
dynamic- response factor increases from 1.15 to 1.27 as the number of 
broken filaments is increased from one to infinity. A somewhat lower 
dynamic- response factor is obtained when a hole is suddenly caused in 
the filament sheet. 


INTRODUCTION 


Structures fabricated from fine filaments that are wound, woven, 
or plied are becoming prevalent in flight applications. Many solid- 
propellant rocket-motor cases, for instance, are being constructed by 
winding re sin- coated glass filaments on a mandrel. The high-strength 
glass filaments carry the pressurization loads and the resin forms a 
matrix which produces a unitized efficient material. Other applica- 
tions make use of the good foldability of coated fabrics to package 
large, low- density structures into small volumes until their erection 
by mechanical means or inflation is desired. 

One of the necessary factors in the rational design of any struc- 
ture is a knowledge of the behavior of stresses in the neighborhood of 
discontinuities such as holes and reinforcements. Whereas much infor- 
mation is available about stress concentration in materials which can 
be considered as continuums, little is known about such stress behavior 
of filamentary configurations . 

Whenever one or more fibers is suddenly broken in a fabric under 
stress, the load in the broken fiber or fibers must be transferred 
through the matrix to the adjacent fibers in order to restore equilib- 
rium. Of interest is not only the resulting static stress state but 
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also the dynamic overshoot which occurs during the transient phase. Both 
of these results are obtained in this paper for several types of "cutouts" 
in the simple case of an infinite flat sheet of parallel filaments 
stressed in uniform tension along the direction of the filaments. 

The model treated is that which is common in shear- lag analyses; 
that is, it is composed of tension- carrying elements connected by purely 
shear- carrying material. The static problem is solved first and the 
details of the dynamic analysis are relegated to the appendix. The 
results are essentially exact within the framework of small- deflection 
elasticity theory. 


SYMBOLS 

a,b major and minor axes of ellipse 

d filament spacing 

EA extensional stiffness of a filament 

Oh shear stiffness of the matrix 

K c stress- concentration factor for a circular hole 

K e stress- concentration factor for an elliptical hole 

Kj, stress- concentration factor for r broken filaments 

1^ load in nth filament for influence- function solution 

m mass per unit length associated with a filament 

n,m indexes 

p applied force on each filament at infinity 

p_ load in nth filament 

n 

P n dimensionless load in nth filament 

r number of broken filaments 

R( ) 


real part of variable 


ro o ui h t~» 
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u 


n 


U 


U 


n 


'n 


x 

y 

z 

^r 

0 

I 

o 

T 


Laplace transform variable 
time 

displacement of nth filament 

dimensionless displacement 
dimensionless displacement of nth filament 

displacement of nth filament for influence- function solution 
coordinate parallel to filaments 
coordinate normal to filaments 
complex variable 

dynamic- response factor for a circular hole 
dynamic- response factor for r broken filaments 
transform variable 

dimensionless coordinate parallel to filaments 
dimensionless stress 
dimensionless time 
elliptical coordinates 


yg _ *( ) , 9( ) 

dx^ <$y^ 

an asterisk denotes Laplace transform in time 


ANALYSIS 


... ® ie conf,i guration under consideration is shown in figure 1 together 
wxth the coordinate and notation systems. The filaments are separated 
by a constant distance and are numbered from -oo to oo from the bottom 
upward. The coordinate along the filaments is denoted by x and the 
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displacement of the nth filament at location x and time t is given 
by u^x, t). Similarly, the force in the nth filament (positive in ten- 
sion) is called p n (x,t) and is given in terms of u n by 


3u n 


P = EA 
n 9x 


( 1 ) 


where EA is the extensional stiffness of the filament. The shear 
force per unit length in the bay between the nth and (n + l)st fila- 
ment is Gh(u n+ 3 _ - %)/d where Gh is the shear stiffness of the 

matrix. Equilibrium of an element of the nth filament then requires 


EA 


dx 2 


+ 


a . (^+1 " 2u n + ^-i) 



( 2 ) 


where the assumption has been made that the mass per unit length m 
associated with the nth filament is concentrated at that filament. 


In figure 1, filaments 0 and 1 are shown broken at x = 0 and the 
remainder are shown intact. In general, for r broken filaments, let 
0 = n = r - 1 denote the broken filaments. The appropriate boundary 
conditions are: 


p n (°,t) = 0 
u n (0, t) = 0 


(0 = n = r - 1) 

> 

(n < 0 or n ^ r) 

) 


( 3 ) 


For x large, of course, the force in each filament approaches the 
uniform applied force which is denoted as p. Thus 


P n (+o°,t) = P 


(4) 


For the time- dependent problem, the following initial conditions 
are required: 


p n (x,0) = p 


Su n 

C)t 


(x,0) 





J 


( 5 ) 
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Nondimens ionalizat ion 

In order to obtain a convenient form for the problem, let 


Pn = P p : 


n 


^ = P 


d 


EAGh Un 


x = 


EAd 

Gh 


. ,/md 

t = T 

(Gh 


( 6 ) 




From these equations the following partial- differential — difference equa- 
tion is obtained: 


a 2 u. 


n 


+ U , i - 2tT + U, 


8% 


n 


a*' 


with boundary conditions 


n+1 


n T n-1 - ~2 
or 


(T) 




U n (0,T) =0 (n < 0 or n r) 

P n (0, t) =0 (0 “S n S r - 1) } (8) 

P n ( ±00 -’ T ) = 1 


and initial conditions 


p n U,o) = i 


du. 


- (1,0) = 0 


J 


( 9 ) 


The dimensionless forces and displacements are related by 

8u n 

P n^' T) =ST 


( 10 ) 
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Solution of Static Problem 

The boundary- value problem for static loading is constituted by w 

equation (7) with the right-hand side set equal to zero and boundary 
conditions (8). The solution is complicated by the fact that the bound- 
ary conditions at x = 0 are mixed; that is, they apply to neither U n 
nor P n solely. The following approach is convenient for overcoming 
this difficulty. 

Influence- function technique .- Consider a filament sheet which has 
no applied edge load and in which all filaments but the zeroeth one are 
intact. Displace the end of the zeroeth filament a unit amount, main- 
tain zero displacement at x = 0 of all the other filaments, and denote 
the resulting forces and displacements by Ln(£) and V n (£)* This set 

of influence functions can then be superposed to obtain the actual prob- 
lem in the following manner: 


P n(^) = 1 + \ L n -mU) U m (0) 

L \ 

m= - °° 




U n (D - I + 2^ V n-mU) U m (0) 

m= - 00 

But boundary conditions (8) yield first 

r-1 

P n^) ~ ^ + ^n-m^ ) ^m^^ 

m=0 
r-1 


u n( 0 =l+^ V n _ m (|) U m ( 0 ) 


m=0 


J 


since U m (0) = 0 for other values of m, and second 


un 


r-1 - 

0 = 1 + V Ln-m(O) U m (0) (0 < n < r - l) (12) 

L 1 

m=0 


which is the specification of the boundary conditions on the loads. 


IY) O \j\ H tr' 
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Equations (12) constitute a set of r equations for the r unknowns 
U m (0). They can "be solved and substituted "back in equations (ll) to yield 

the entire solution. First, however, the L n values must be determined. 

Determination of the influence functions Lq.- The problem can be 
stated for | > 0 as ~~ “ 


d 2 V. 


n 


d|' 


+ v n+l * 2v n + v n-l = 0 


(13) 


with the conditions 


v n (0) = 1 

V„(0) = 0 


dVi 


- (°°) = o 




In order to solve this problem, let 


v(s,e) = 


V 


(n = 0) 
(n/0) j> 


v n (D 


-in0 


n= - to 


(HO 


(15) 


or, inversely. 


v n (i)=ij^ v(i,e) e ine de 


Then, multiplying equation (13) by e'* 110 


(16) 


a 2 ? 


at 


2 - k sin 2 | V = 0 


and summing over all n gives 

(17) 


Similar treatment of the boundary conditions (l4) yields 

v(o,e) = 1 

^ (00,9) = 0 


(18) 
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The solution satisfying equations (17) and (l8) is 

-2 


V = e 


sing 


Thus, from equation (l6) 

p it 

v n (l) = i J COS ne 


-2|sin— 
e 2 d6 


(19) 


which can he expressed in terms of Bessel and Weber functions of imagi- 
nary argument. For the present purposes, however, the reduction is not 

dV n 

as 


necessary since attention will be centered on 1^(0) 


(0) which 


is given simply by 


w 


L 


L n (°) 


(ito 2 - l)it 


( 20 ) 


Stress- concentration factors .- Inspection of the problem indicates 
that the maximum force occurs at | = 0 in the first intact filament 
adjacent to the broken ones. Thus, the stress- concentration factor K r 
for r broken filaments is given by P r (0) which equals P_q(0) by 
virtue of symmetry. Now, 


r-1 


P r (0) 


= 1 + 


Vm(°) 


Um(0) 


m=0 


( 21 ) 


Solution of equations (12) for the U m (0) and substitution into 

equation (2l) yields the stress- concentration factor. This process has 
been carried out for an r of 1 to 6 and the results are given in the 
following table 


r 

K r 

1 

V5 

2 

8/5 

3 

6^/35 

if 

128/63 

5 

512/231 

6 

l, 02 if/if 29 


OUl H 
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Inspection of these values show that they can be written as 


- it ’ 6 • 8 • I • ( 2r + 2 ) 

^ 3 . 5 . 7 . . . (2r + 1) 


( 22 ) 


Although this result has not been established in general, its correct- 
ness for the first six values lends credence to its validity for all 
values of r. 


Solution of Dynamic Problem 

If a Laplace transform is taken of the time- dependent differential 
equation (eq. (7)) and boundary conditions (eq. (8)), the resulting equa- 
tions are similar in form to the static equations and the same type of 
approach can be used for their solution. The details of the solutions 
are considerably more complicated and are contained in appendix A. The 
resulting timewise variation of the stress- concentration factor is shown 
in figure 2 for 1, 2, and 3 broken filaments . Solutions for greater 
numbers of broken filaments were not obtained because of the increasing 
difficulty of calculation and because of the existence of an apparent 
upper limit on the dynamic overshoot. This upper limit is discussed in 
the next section. 


DISCUSSION OF RESULTS 


As can be seen from figure 2, the stress- concentration factor 
exhibits an oscillation that decays in a few cycles to the steady- state 
value, as the energy is carried away to infinity by wave motion. The 
oscillations for r = 1 and r = 2 are fairly simple and are similar 
to those of a one- degree- of- freedom system. For r = 3, the oscillation 
is more complex, an apparent second-mode component appearing in the time 
history. In all cases, the first peak is the largest onej the value of 
the stress at this peak determines the dynamic overshoot and is the 
principal result to be extracted from the dynamic analysis. 


Ifynamic- Response Factor 

The dynamic-r response factor q r is defined as the ratio between 

the maximum stress and the static stress. Values for 1, 2, and 3 broken 
filaments are given in the following table: 



r 

*l r 

1 

1.15 

2 

1.19 

3 

1.20 


The dynamic overshoot thus apparently increases with increasing number 
of broken filaments. The overshoot can therefore he reasonably expected 
to be the highest in the limit as the number of broken filaments approaches 
infinity. 

Results for an infinite number of broken filaments .- As has been 
noted, the greater the number of broken filaments, the greater the dif- 
ficulty of solution; but the limiting case itself is readily amenable 
to analysis and is treated in appendix B. 

The analysis in appendix B deals with the so-called continuous 
stringer sheet which is an orthotropic medium with finite extensional 
stiffness in the longitudinal direction, infinite extensional stiffness 
in the transverse direction, and finite shear stiffness. Its behavior 
is governed by the nondimen si onali zed differential equation 

V^U = — (23) 

£>r 2 

which is obtained either by direct derivation or by replacing the second- 
order difference in equation ( 7 ) with its appropriate derivative 
equivalent . 

The breaking of an infinite number of filaments is accomplished by 
placing a finite- length slit in the sheet. As is well-known, this pro- 
cedure leads to an infinite stress- concentration factor in the static 
problem, the stress varying as the inverse one-half power of the dis- 
tance from the end of the slit. The stress behaves similarly in the 
dynamic problem. Thus, the dynamic behavior can be studied by finding 
the variation of the strength or magnitude of the stress singularity 
with time. This central result of appendix B is shown in figure 3- As 
can be seen, the strength C exhibits discontinuities in slope. These 
discontinuities arise from the reflection of waves issuing from the ends 
of the slit. Only the first two regimes are shown in figure 3 together 
with the eventual static value. Results for subsequent regimes are very 
difficult to obtain and apparently would not contribute any greater 
stresses. The resulting dynamic- response factor is 

t] = 1-27 
00 

which should be an upper limit on T] r . 
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Hole in a stringer sheet .- In the foregoing sections, only slits 
have been treated. Of interest also is the case in which a hole is 
punched out of the material. Because of the usefulness of conformal 
mapping techniques, the static problem is easily analyzed for a large 
variety of hole shapes for a stringer sheet. In particular, for ellip- 
tical shapes, the stress- concentration factor is derived in appendix C 
to be 

K e = 1 + e ( 2 k) 

where e is the ratio between the transverse and longitudinal dimen- 
sions of the elliptical hole in the nondimens ional coordinate system 
pertaining to equation (25). For example, a stringer sheet with an 
elliptical hole that transforms into a circle in the nondimens ional 
coordinate system has a static stress- concentration factor of 2 . 

In this latter case of a circular boundary, the dynamic problem is 
also tractable. The analysis is given in appendix C and the variation 
of K c with time is shown in figure The resulting dynamic- response 
factor is 

q c = 1.08 

which is considerably less than that for the slit. 


CONCLUDING REMARKS 


The results derived in this paper lead to the apparent conclusion 
that the highest dynamic- response factor applicable to the stresses 
acting in the neighborhood of suddenly induced discontinuities is I.27. 
This value was obtained for the case of a long slit in a sheet of closely 
packed filaments. Either reducing the number of broken filaments or 
opening the slit into a hole can be expected to decrease the dynamic- 
response factor. 

The analysis was based on elastic, small- deflection theory of a 
two-dimensional medium. In actuality, filamentary sheets usually have 
a large number of filaments through the thickness and, whereas the fila- 
ments themselves may follow Hooke’s law very well, the matrix or the 
weave introduces large nonlinear effects in the form of plastic defor- 
mations, large deflections, and fiber straightening. These factors 
should act to reduce the stress concentrations and, at least, not 
increase the dynamic overshoot. In view of the smallness of the dynamic 
overshoot relative to the possible effects of these other factors, future 
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theoretical work on this subject would seem to be best devoted to ana- 
lyzing better models of various types of filamentary construction. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va., March 20 , 1961. 
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APPENDIX A 


DYNAMIC BEHAVIOR OF SHEET WITH SUDDENLY BROKEN FILAMENTS 


Let the Laplace transform in time of U n (|,r) he denoted by 
U n (|,s) where s is the transform variable. Then, taking transforms 
of equations ( 7 ) and (8) yields 

a 2 u* 




- + u n+l - (2 + s2 ) u n + u n-l = 


(Al) 


U n (0,s) = 0 
P*(0,s) = 0 

/ .-j 

sf -i 




(n<0 or n ^ r) 

(0 = “ = r • 1) l (A2) 


J 


where in initial conditions ( 9 ) the condition on P n has been converted 
to the condition U n (£,0) = £. 

Again, use can be made of the unit solution to write the solution 
for transformed loads and displacements in the form 

r-1 


Pn(S»a) = | + Y, L °-“ (l ’ s) ^ (0 ^ s) 


U n (l> s) — "g + \ ^ V n _ m (|,s) s) 


m=0 

r-1 


m=0 

from which the following equations are obtained 

r-1 

1 


(A3) 


J 


0 = - + L*_ m (0,s) U*(0,s) (0 = n = r - l) (Ah) 


m=0 



14 


and from which the unknown 1^(0, s) can be determined. Also, the trans- 
formed load in the first intact filament is 


r-1 

1 


P r (0, s) - g" + \ L r _ m (0, s) U m (0,s) 


s ’ ) r-m' 

m=0 


(A5) 


The transform of the stress- concentration factor for r broken 

filaments is set equal to this transformed load. 

ay* 

As before, Lj(|,s) = — (i,s) where vj(|,s) satisfies 

d| 


2 ,,* 




* 1 ‘ 


a + Ci - (2 + s Z K + Cl = ° 


and the conditions 


V*(0,s) = 1 

'S' 

It 

C3 

V*(0,s) = 0 

^ (»,s) - 0 

(n ^ 0) 


j 


In order to solve these equations, let 

V"( 6,s,e) - \ V*(|,s) e- ln0 


n=- oo 


or, inversely. 


(A6) 


(A7) 


(A8) 


Vn(l,s)=i / V*(|,s,0) e in0 d9 

2lt J -% 


(A9) 


,-inG 


Then equations (A 6) and (A7) become, after multiplying by 
summing. 


and 



ro o ui m r 
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a 2 ?* 

as 2 


- ^4 sin 2 | + s 2 )v* = 0 


ai 

The solution for | > 0 is 


v*(o,s,e) = 1 

= 0 




-r& 

V = e 


-ft. sin 2 i 


+ s 2 I 


which yields 


V*(|,s) = i y cos n6 e * 


-|/4sin 2 -l+s 2 | 


d9 


or, finally 


1^(0, s) ^ ~ J' cos n0 ^ sin 2 -| + s 2 d0 


(A10) 


(All) 


(A12) 


(A13) 


Solving equations (A4) and substituting into equation (A5) gives 
the transform of the stress- concent rat ion factor. For 1, 2, and 3 broken 
stringers, this procedure yields 


K-.(s) 


-t* t* n 
L 0 ~ L 1 

v 4 ' 


4(b) -i 


I-S - L 2 
* * 
L 0 + L 1 


K^(s) 


(L$ - Lj)(L$ - 1$) - (ij - L %f 
L$ 2 + L^L2 - 2L l 2 


(Alta) 


(Al4b) 


(Alta) 


where, for brevity, the functional dependence indicated in equation (A13) 
is omitted. 
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The task remains to evaluate L^(0, s) as a function of s and to 

take the inverse transform of each of equations (Al4) . This inversion 
requires integration in the complex s-planej a study of the behavior of 
L* ^and, therefore, of K* j as functions of the complex variable s is 

necessary. Such a study shows that has branch points at s — I2i_,0. 

(Actually, L-q can be written in terms of complete elliptic integrals 
with modulus 



These elliptic integrals have branch points at k = ±1. The form of the 
modulus also guided the choice of conformal mapping used in the sequel.) 
But the square root in equation (A13) can be thought of as behaving 
essentially as s for large values of s. Therefore, the branch cuts 
need not extend to infinity and l£ can be made single valued by placing 

a branch cut along the imaginary axis between -2i and 2i. 


The inversion integral is 


Mt) 


1 

2rti 


p 7+ieo 
y— i 00 


K*(s)e ST ds 


Since the integrand satisfies Jordan’s lemma and the denominators of 
K* have no zeros except possibly on the branch cut, the integral path 

can be closed around infinity on the left-hand side and shrunk to the 
contour C around the branch cut. Thus, 


K r (r) 



(A15) 


A closed- form evaluation of this integral is probably impossible 
so that a series evaluation is desirable. One method which is very 
suitable is as follows: 


Let 


s 


= z - 


1 

z 


(ai6) 


This maps the entire s- plane outside the branch cut into the interior 
of the unit circle. The counterclockwise contour around the branch 
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cut maps into a clockwise contour just inside the unit circle. 
Furthermore, 


^(°, z ) = - i J'. 

J 0 


n fi 

cos n9 ' 


p II 

2z^cos 9 + 


d0 


which is expandable in a power series in z that is convergent inside 
the unit circle. Thus, 


l£(0,z) 


1 

z 


00 


E 

k=0 



(A17) 


where 



is the binomial number. 

Making the change of variable gives 



(ai8) 


The integrand has no singularities within the path of integration except 
at the origin which i# an essential singularity. The value of K r (r) 

is therefore given simply by the residue of the integrand at z = 0. 

The determination of this residue involves finding the coefficient of 
the zeroeth power of z in the expansion of 



z^ + 1 
z 2 - 1 



T 


e 


The series for the bracketed part can be determined from equations (A17) 
and (Al4) . It involves only positive even powers of z. So only the 
negative even poVers in the expansion of the other part must be sought. 
If 



n = -co 


(A19) 
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there results 


C 0 - -l 


C_ n - Jq(2t) + 2J2(2t) + ♦ • • + 2J n _2(2r) + J n (2r) - 1 (n = 2,4,6, . . .) 


(A20) 


where J n is the Bessel function of the first kind. 

Proper manipulation and evaluation of the series gives the results 
plotted in figure 2. Good convergence is obtained and the results are 
quite accurate. 


^ H IAO C\J 
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APPENDIX B 

DYNAMIC BEHAVIOR OF A SUDDENLY SLIPPED STRINGER SHEET 



The required stress distribution can be written 

du 


a = 1 + 


dx 


(Bl) 


where u satisfies the following dimensionless boundary- value problem 


V 2 u - 

£>t 2 


(B2) 


u(x,y,0) = (x,y,0) = 0 




Also, u is 
running from 


^ (o,y,T) = -l 


(-1 < y < 1) [ 


(B3) 


J 


regular everywhere in the x,y plane except on the slit 
y=-l to y=l and approaches zero at infinity. 


Rather than solve this boundary- value problem directly, the use of 
an analogy to obtain the stress behavior near the ends of the slit is more 
convenient. The problem as stated is in exactly the same mathematical 
form as the problem of determining the perturbation velocity potential 

in steady supersonic flow (with Mach number of (2) over a thin rectan- 
gular wing (of span 2) with a very long chord. The time coordinate r 
corresponds to the coordinate in the stream direction and the x and 
y coordinates correspond to the coordinates in the crossflow plane. 

Zero time is equivalent to the leading edge of the wing which has a 
constant slope of -1 in the stream direction. Finally, the stresses 
along the y-axis are the same as the upwash velocities in the tip 
regions in the plane of the wing. 


The existence of this analogy allows the use of the well developed 
methods of supersonic wing theory; that of Eward in reference 1 is 
particularly applicable. This method obtains the upwash velocities in 
terms of explicit quadratures which are expressible in terms of tabu- 
lated functions for the regions close to the leading edge. 







20 


The actual detailed analysis by Eward's method is somewhat tedious 
and only the resulting formulas of interest are given here. The analysis 
shows that the stress at the y-axis in the neighborhood of the tip y = 1 
is of the form 


a = + 0 ( 1 ) 

\ly - 1 

where C(t) is the strength of the singularity and is given by 


0( T ) = - |/t 

Jt 


C(t) = f~T~2 (B* - K') + ft 


$ £ + K'Z(v) 
2 K 


(0 < t < 2) 


(2 < r < b) 




(B*0 


> (B5) 


where K and E are complete elliptic integrals of the first and sec- 
ond kind with modulus 


I 

1 


k 


r + 2 


and the primes refer to functions of the comodulus 


k’ = 


t - 2 


t + 2 

The quantity v is determined from the equation 


sn v = 


t + 2 
2t 


where sn v is Jacobi’s elliptic sine. Finally, Z(v) is Jacobi's 
zeta function. (See ref. 2.) 

The results are plotted in figure 3* As can be seen, the maximum 
value of C is reached at t = 2. The asymptotic value for large t 
is obtained from the solution of the static problem: 

= 0 

(0,y,r) = -1 (-1 < y < l) 


m n vji 



ro o vji h f 
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The solution can he seen to he 


u(x,y) = R l |/(x + iy) 2 + 1 - x 


whence the stress is 


or = R 


x + iy 


y(x + iy)' 


+ 1 , 


In particular, the stress at the y-axis is 

a(0,y) = 0 


cr(0,y) = 




^ - 1 


( |y| < i) 
(y > i) 


C(co) 


f2 


so that 
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APPENDIX C 


STRESS CONCENTRATIONS IN A STRINGER SHEET WITH A HOLE 


The stress in an infinite stringer sheet with uniform stress in the 
x- direction at infinity and with a suddenly induced hole can. be written 


a - 


du 

dx 


(Cl) 


where u satisfies 


and the following initial conditions 

u(x^y,0) = x 

du 


(C2) 




C^T 


{■x.,y,0) = 0 


(C3) 


and boundary conditions 
at the boundary of the hole: 


— = 0 

dn 


(C4a) 


and at infinity: 


du 

dx 


->1 


(C4b) 


These equations are solved in the sequel for two cases: the static 
stresses around an elliptical hole and the dynamic stresses around a cir- 
cular hole. 


Static Problem 

For the first case, let 


x + iy = c cosh(cp + i^) 


ro O VJI H M 
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Then the ellipse 



maps into 


cp = cp 0 


where 


c cosh cpQ = a 
c sinh cpQ = b 

The differential equation and boundary conditions 

+ ^£u = 0 
dcp 2 dip 2 

O'*) =° 

u(cp,\|/) ~ -2 e^cos i|r 


become 


(<P ^ <Pq) 0 < t < 2 jt) 


(0 < t < 2 t r) 

(cp » cp 0 ; 0 < \|r < 2*) 


This problem has the solution 


u(cp,\Jr) = — 2 — (a cosh cp - b sinh cp) cos \|/ 
a - b 

Therefore, 

a = 1 / b cosh cp sinh cp \ 

a ~ ^ \ cosh 2 cp - cos 2 \j/ / 

The maximum stress occurs at \|r = cp = cp Q (at the boundary of the 
hole on the y-axis) and is 



As derived, this result is valid only for a > b. But a similar deriva- 
tion for a ^ b yields the same result. Therefore, the stress concen- 
tration factor is 


Kg = 1 + e 


(C5) 


where e is the ratio of the transverse and longitudinal axes of the 
ellipse . 



2k 


Dynamic Problem 


For the dynamic problem with a circular hole of unit radius, let 

(c6) 

t 

y = r sin 0 


x = r cos 0 


and take the Laplace transform in time. (See appendix A.) The following 
boundary- value problem results 


H* + i + J_ HI - sV = -sr cos 0 
Sr 2 r SF r 2 Se 2 


(C7) 


The solution is 


Su* 

Sr 




(1,0, s) = 0 


u*(r, 0, s) 


r cos 0 


(r » l) 


j 


u*(r, 0,s) = 


COS 0 


r - — 


i %( sr ) 


s h( s ) 


(C8) 


(C9) 


where is the modified Bessel function of the second kind and the 

prime denotes differentiation. 

The transform of the stress at x = 0, y = 1 is 


* 


a 



1 1 Kl(a) 

S s 2 Ki(s) 


(CIO) 


The stress-concentration factor is given by the inverse transform of 
equation (CIO). Thus, 


K c (t) = 1 - 


2jti 


r 7+1 “ K l(s) e ST ds 
7-ic 


Ki(s) 


(Cll) 


The integrand has a branch point at s = 0 and can be made single 
valued by placing a branch cut along the negative real axis. The 


5 

c 


1 — I I — I 
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0 

p 


integrand also das two poles located symmetrically in the left hand half- 
plane which arise from the zeros of Kq. By taking account of the resi- 
dues at these poles, the integration path in equation (Cll) can be warped 
into the path around the branch cut. This latter integral can then be 
rewritten in real form to yield after some manipulation 


SqT S 0 T 

K c (t) = 2 - -2 - e --_ --g 

1 + Sq 1 + s 0 


+ 



e -Ta) dm 



(C12) 


where s 0 and sq are the (conjugate) zeros of Kq(s), and Iq and Ip 
are modified Bessel functions of the first kind. 


The value of sq was found approximately by interpolation by using 
the zeros of Kq/^j ^5/2* approximate value was then 

refined by means of the Newton iteration method by using the series 
expansion for Kq. The result is 


s 0 = -0.64355 + 0.501l8i 


The integral in equation (C12) was evaluated numerically after 
making the substitution 


v 


e 


-2m 


Ten intervals in the range 0 < v < 1 were used together with Simpson’s 
rule . 


The results are accurate except for the smallest values of t. The 
stress- concentration ratios for these small values were obtained by per- 
forming a term-by-term inversion of the asymptotic- series expansion of 
the right-hand side of equation (CIO). The resulting small-time approxi- 
mation is 


K c (t) 


= 1 + T - 


4 



( T « 1) (C13) 


The results are plotted in figure 4. 
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